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Abstract 

We study the b —* S7 decay including the next to leading QCD corrections in the 
two Higgs doublet model with flavor changing neutral currents at the tree level. We find 
the constraints to the flavor changing parameters of the model, using the experimental 
results on the branching ratio of B — » X s j decay, provided by the CLEO Collaboration, 
the restrictions coming from the AF = 2 (F = K,D,B) mixing and the p parameter. 
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1 Introduction 



Rare B meson decays are one of the most promising research area in particle physics and lie 
on the focus of theoretical and experimental physicists. In the Standard Model (SM), they are 
induced by flavor changing neutral currents (FCNC) at loop level and therefore sensitive to the 
gauge structure of the theory. From the experimental point of view, they play an outstanding 
role in the precise determination of the fundamental parameters of the SM, such as Cabbibo- 
Kobayashi-Maskawa (CKM) matrix elements, leptonic decay constants, etc. Furthermore, these 
decays provide a sensitive test to the new physics beyond the SM, such as two Higgs Doublet 
model (2HDM), Minimal Supersymmetric extension of the SM (MSSM) etc. Among the 
rare B decays, b — > S7 has received considerable interest since the branching ratios (Br) of the 
inclusive B — > || and exclusive B K*j H have been already measured experimentally. 
Recently, the new experimental results for the inclusive b — > 57 decay are announced by CLEO 
and ALEPH Collaborations ||. Therefore, the b — > s'j decay is under an extensive investigation 
in the framework of various extensions of the SM, in order to get information about the model 
parameters or improve the existing restrictions. 

It is well known that the FCNC are forbidden at the tree level in the SM. This restriction 
is achieved in the extended model with the additional conditions. 2HDM is one of the simplest 
extensions of the SM, obtained by the addition of a new scalar SU(2) doublet. The Yukawa 
lagrangian causes that the model possesses phenemologicaly dangerous FCNC's at the tree 
level. To protect the model from such terms, the ad hoc discrete symmetry || on the 2HDM 
scalar potential and the Yukawa interaction is proposed and there appear two different versions 
of the 2HDM depending on whether up and down quarks couple to the same or different scalar 
doublets. In model I, the up and down quarks get mass via vacuum expectation value (v.e.v.) 
of only one Higgs field. In model II, which coincides with the MSSM in the Higgs sector, the up 
and down quarks get mass via v.e.v. of the Higgs fields Hi and if 2 respectively where Hi(H2) 
corresponds to first (second) Higgs doublet of 2HDM 0. In the absence of the mentioned 
discrete symmetry, FCNC appears at the tree level and this model is called as model III in 
current literature |7], |8], [| . A comprehensive phenemological analysis of the model III was done 
in series of works |7], n|. In particular, from a purely phenomenological point of view, low 
energy experiments involving K° — K°, B° — B°, Kl — > etc, place strong constraints on 
the existence of tree level flavor changing (FC) transitions, existing in the model III. 

In the present work, we examine the b —>■ S7 decay in the model III, taking the next to 
leading (NLO) QCD corrections into account, in a more detailed analysis compared to one 
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existing in literature (see |7], |10|). Further, we obtain the constraints for the neutral couplings 



Cmu £.Nbb an d £,Ntc with the assumption that £jv cc , £,Nsbi £,n ss an d the other couplings which 
include the first generation indices are negligible compared to former ones (for the definition of 
£,N,ij see section 2). Our predictions are based on the CLEO measurement B — > X s ^ and the 
restrictions coming from the AF = 2 (F = K,D,B) mixing and the p parameter \W\. Note 
that NLO QCD corrections to the b — > 37 decay in 2HDM (for model I and II) were calculated 

DJ Pi- 



rn 



The paper is organized as follows: In Section 2, we present the NLO QCD corrected Hamil- 
tonian responsible for the b — > s'j decay in the model III and discuss the effects of the additional 
Left — Right flipped operators to the decay rate. Section 3 is devoted to the constraint analysis, 
more precisely to the the ratios 

^Nbb 



and our conclusions. 



2 Next to leading improved short-distance contributions 
in the model III for the decay b —> 57 

Before presenting the NLO QCD corrections to the b — > s'j decay amplitude in the 2HDM 
(model III), we would like to remind briefly the main features of the 2HDM. The Yukawa 
interaction for the general case is 



£y = vfjQiL^lUjR + VuQiLfaDjR + &QiL<hUoR + Zf-QiLfaDjR + h.c. 



(1) 



where L and R denote chiral projections L(R) = 1/2(1 =p 75), fa for i — 1,2, are the two scalar 
doublets, r]fj D and are the matrices of the Yukawa couplings. For convenience we choose 
fa and 02 in the following basis: 



V2 





v + H° 



+ 



( V2x 
{ ^X 



J_( V2H+ 



(2) 



where the vacuum expectation values are, 

<^>=^= 2 {l) ;<0 2 >=o. 

This choice permits us to write the FC part of the interaction as 

£y,fc = £ijQiLfaU jR + ^Q iL faD jR + h.c. 



(3) 



(4) 



with the following advantages: 
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• doublet 0i corresponds to the scalar doublet of the SM and H Q to the SM Higgs field. This 
part of the Yukawa Lagrangian is responsible for the generation of the fermion masses 
with the couplings r) u,D . 

• all new scalar fields belong to the <p2 scalar doublet. 

The couplings £ U,D are the open window for the tree level FCNC's and can be expressed for 
the FC charged interactions as 

£ci = ^neutral VcKM , 

idi = VCKM ^neutral , (5) 

where Cneutrai Q * s defined by the expression 

= (v u L > D y^ D v u R D . (6) 

Here the charged couplings appear as a linear combinations of neutral couplings multiplied by 
Vckm matrix elements. This gives an important distinction between model III and model II 

(I)- 

After this preliminary remark, let us discuss the NLO QCD corrections to the h — > 57 
decay in the 2HDM for the general case. The appropriate framework is that of an effective 
theory obtained by integrating out the heavy degrees of freedom, which are, in this context, 
t quark, W , H , Hi, and H 2 bosons, where H ± denote charged ,Hi and H 2 denote neutral 
Higgs bosons. The LLog QCD corrections are done through matching the full theory with the 
effective low energy theory at the high scale /1 = and evaluating the Wilson coefficients 
from mw down to the lower scale fi ~ 0{mb). Note that we choose the higher scale as /1 = mw 
since the evaluation from the scale /1 = m#± to fi = my? gives negligible contribution to the 
Wilson coefficients. Here we assume that the charged Higgs boson is heavy due to theoretical 



analysis of the b — > 57 decay (see [[11], [T 

The effective Hamiltonian relevant for b — > sj decay is 

. G F 



H 



eff 



^V tb V*Y,C^)0^) , (7) 



where the Oi are operators given in eq. (§) and the Cj are Wilson coefficients renormalized at 
the scale \i. The coefficients are calculated perturbatively and expressed as functions of the 
heavy particle masses in the theory. 



1 In all next discussion we denote Neutral as £zv ■ 
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The operator basis depends on the model used and the conventional choice in the case of 
SM, 2HDM model II (I) and MSSM is 

01 = {s Lo r( ll c L p)(c L p'fb La ), 

2 = {s Lo fl ll CL a ){c L p'fb L p), 

q=u,d,s,c,b 

04 = (sLalpbLp) J2 (QLp^QLa), 
q=u,d,s,c,b 

C 5 = (sLal^b La ) J2 (wAfl/?), 
q=u,d,s,c,b 

q=u,d,s,c,b 

g 

07 = j^s a a lxu (m b R + m s L)b a J 7 t 11 ', 

8 = ^s a T« p <j^{m b R + m s L)bpQ a »\ (8) 

where a and (3 are SU (3) colour indices and T^ v and are the field strength tensors of the 
electromagnetic and strong interactions, respectively. 

In our case, however, new operators with different chirality structures can arise since the 
general Yukawa lagrangian includes both L and R chiral interactions. The conventional operator 
set is extended first adding two new operators which are left-right analogues of 0± and O2, 
namely 

09 = (s La 'y li c L i3)(c R i3ry fi b Roi ), 

0io = {sL a l^c La ){c R ^b R p), (9) 
Further we need the second operator set O^—O'^ which are flipped chirality partners of 0±— O w : 



o[ = 


(sR a 'y ll C R f3)(c Rf 3>y%Ra) , 


0> = 




O's = 


q=u,d,s,c,b 


01 = 






q=u,d,s,c,b 


05 = 


(sRa-f^Ra) {QLpl^Lp), 




q=u,d,s,c,b 


0; = 


q=u,d,s,c,b 


0; = 


g 

167r2 s a cv(m 6 L + m s R)b a FL ll/ , 
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0' 8 = jl^s a T^v(m b L + m s R)bpg a ^, 
0' g = (s i?a 7 M c J?/3 )(c L/3 7 M 6 La ) , 

O[ = (s i? , a 7 A1 c fia )(c L/3 7 M 6 L /3) • (10) 
This extended basis is the same as the basis for SU{2) L x SU(2) R x £7(1) extensions of SM 



nj. Note that in the SM, model II (I) 2HDM and the MSSM, the absence of 0' 7 and 0' 8 are a 
consequence of assumption m s /m h ~ 0. 

In the calculations, we take only the charged Higgs contributions into account and neglect 
the effects of neutral Higgs bosons for the reasons given below: The neutral bosons Hq, Hi and 
H 2 are defined in terms of the mass eigenstates H ,/i and A as 

H = (H cosa — hosina) + v , 
Hi = (h cosa + Hasina) , 

H 2 = A, (11) 

where a is the mixing angle and v is proportional to the vacuum expectation value of the 
doublet 0i (eq. @). Here we assume that the massess of neutral Higgs bosons ho and A are 
heavy compared to the b-quark mass. The neutral Higgs scalar ho and pseduscalar Aq give 
contribution only to Cj for b —>■ s'y decay With the choice of a = 0, Cj Q and can be 
expressed at mw level as 



Qi 



i=d,s,b 



i rrii rat, 



cf°(m w ) = {v tb v: s y l E c^,,; (a'., • (12) 



i=d,s,i 



i rrii mb 



where rrii and Qi are the masses and charges of the down quarks (i = d, s, b) respectively. Here 
we used the redefinition 



t u ' D = • (13) 

Eq. (0) shows that neutral Higgs bosons can give a large contribution to CV, which does not 
respect the CLEO and ALEPH data . At this stage we make an assumption that the couplings 
£§ is (i = d, s, b) and £®db are negligible to be able to reach the conditions ^Nbb^Nis << 1 an d 
^Ndb^Nds << 1- These choices permit us to neglect the neutral Higgs effects. 

Now, for the evaluation of Wilson coefficients, we need their initial values with standard 
matching computations. Denoting the Wilson coefficients for the additional charged Higgs 
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contribution with C^{my/), we have the initial values of the Wilson coefficients for the first set 
of operators (eqs.(H), (§)) 

rn t v ts v tb 

1 - V* - Vi 

m t V ts V tb 

+ — (Z u n,u + + 8U^)G 2 (v) ■ (14) 

m t m b ' V t * s V tb 

The explicit forms of the Wilson coefficients in the SM (Cf M (m\y)) is presented in the literature 
T5fl . For the primed Wison coefficients we get, 



C 1 SM G ! 9 :W ('m w ) - , 

(15) 



Cg.sjwimw) = 

C' 7 H (m w ) = — 2 (^ bs ^ + ^J(^ bb + ^ sb ^)F 1 (y) 

m t v ts V tb 



mt m b V ts V t b 

C' 8 H (m w ) = — (C^bsTFj + £n,ss) (^N,bb + €N,sb7T-)Gl(y) 7 
mf V t * s Vtb 

+ -V (Hbs^ + S,J + S, te ^)G a (y) , (16) 
mt m b v ts v t b 

where x = mf/m^r and y = m\jm 2 H ±. The functions Fi(y), F 2 (y), G\(y) and G 2 {y) are given 



as 



7/(7- 5y -% 2 ) i/ 2 (3y-2) 
= 72(y-l)3 + l2^W lny ^ 

y(5y-3) , y(-3y + 2) , 
^ = I%^T7 + 6(y-l)3 

y(-y 2 + 5y + 2) -y 2 
= 2% -1)3 + 45Pl? ,ny ' 

= 4(^iy + %^iy^- (17) 
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In calculations we neglect the small contributions of the internal u and c quarks compared to 
one due to the internal t quark. 

For the initial values of the mentioned Wilson coefficients in the model III (eqs. (0), (|15|) 
and (|I6|)), we have 

cP 2HDM (m w ) = cP SM (m w ) + cP H (m w ). (18) 

Using these initial values, we can calculate the coefficients Cf HDM and Q 12 at any lower 
scale with five quark effective theory where large logarithims can be summed using the renor- 
malization group. Since the strong interactions preserve chirality, the operators in eqs. (|^. |^) 
can not mix with their chirality flipped counterparts eq. fllCf ) and the anomalous dimension 
matrices of two separate set of operators are the same and do not overlap. With the above 
choosen initial values of Wilson coefficients, their evaluations are similar to the SM case. 

For completeness, note that, the operators 5 ,0 6 ,Og and O w (0' 5 ,0'q, 0' g and 0' 1Q ) give 
contributions to the matrix element of b —>■ 57 and in the NDR scheme which we use here, the 
effective magnetic moment type Wilson coefficients are redefined as 

CpM = Cr DM {») + Q d {Cl HDM {») + N c Cl HDM {^)) , 

+ Qu Of* DM {n) + N C ^ C$ HDM (n)) , 
m b m b 

C' 7 eff M = C? hdm {li) + Q d (C' 2HDM (fi) + N c C'i HDM ^)) 

+ Qui— C M V) + N c ^ C> 2HDM (»)) , (19) 

TTlf, 171b 

where N c is the color factor and Q u (Qd) is the charge of up (down) quarks. There is still 
another mixing in the operator set 7 , 0$, Og, O10 (0' 7 , 0' 8 , 0' 9 , 0' 1Q ) Jl4| and we do not take 
into account since the initial values of the Wilson coefficients C\q and C' w are zero in our case. 
The NLO corrected coefficients C 7 HDM (jj) and C 7 HDM (fi) are given as 

Aix 

C' 2HDM U) = C 7 LO > 2HDA \ri + —C' 7 {1)2HDM (v) . (20) 

Atx 

Here 77 = a s (mw)/c( s (fi), hi and a« are the numbers which appear during the evaluation [16]. 
The functions C 7 ° ,2HDM (fi) |L7| and C' 7 L ,2HDM (f/,)axe the leading order QCD corrected Wilson 
coefficients: 

C' 7 LO ' 2HDM (fi) = r ? 16 / 23 C'f™ f (m w ) + (8/3)(r ? 14 / 23 -r / 16 / 23 )C'^ DA/ (m w ) (21) 
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and C>? HDM (/i) is the a s correction to the leading order result that its explicit form can 
be found in [0, pf. c'j 1)2HDM '(//) can be obtained by replacing the Wilson coefficients in 
C\ r 2HDM (ii) with their primed counterparts. The NLO corrected coefficients C§ HDM (fi), 
Cq (/i) and C^ HDM (/j,), Cq HDM (fi) are numerically small at m;, scale therefore we neglect 
them in our calculations. 

Finally, the NLO QCD corrected b — > 57 decay rate in model III is obtained as 

r(6 - 87) = ^r« em |^^| 2 (|C 7 e// (m & )| 2 + |^ e// (m b )| 2 ) , (22) 



32tH 



where a em is the fine structure constant, and m b is b-quark mass. \C% (m b )\ 2 is given in |l I 



\C e 7 ff (m b )\ 2 = \D\ 2 + A+%\C?°' 2HDM (m b )\ 2 

mi 

^c P r> r ( y^tLO,2HDM / \\* ( nLO,2HDM i \ 1 ^<LO,2HDM / \\ i , 

+ -^-.Rel^Cy (m 6 )J ^C/ 2 (m 6 ) - -Cj (23) 



'6 

The functions D and A are [|11|] 



16 ~LO,2HDM, 

'nth 



D = C 7 i0 ' 2MM (m 6 ) + g(c^ 
A = £ i?e{^- 2 ™K)fcf' 2 ™K)) , / y }. (24) 



The explicite expressions for /y, r i; 5^ p and 5^ can be found in (TT|. At this point we would 



like to note that the expressions for unprimed Wilson coefficients in our case can be obtained 
from the results in JTT] by the following replacements: 

W t K ts Vife 

1 - ^* - - Vts 

m t m b V t * s ' Vtb 

To obtain |C 7 e ^(m fe ) | 2 , it is enough to use the primed Wilson coefficients at m w level 
(eq. [16]) since the evaluation of C 7 6 ^(/i) from fi = m w to /i = m b is the same as that ofC 7 e// (At). 
Note that, for model II (model I) Y and XY are 

Y = 1 /tan/3 (l/tan(3) , 
XY = 1 (-l/tan 2 p) . (26) 

3 Constraint analysis 

Now let us turn our attention to the constraint analysis. Restrictions to the free parameters, 
namely, the masses of the charged and neutral Higgs bosons and the ratio of the v.e.v. of the 
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two Higgs fields, denoted by tan (3 in the framework of model I and II, have been predicted in 
series of works |18| . Recently the constraint which connects masses of the charged Higgs bosons, 
rriH± and tanf3, is obtained by using the QCD corrected values in the LLO approximation and 
it is shown that the constraint region is sensitive to the renormalization scale, p fl~3|| . 

Usually, the stronger restrictions to the new couplings are obtained from the analysis of the 
AF = 2 (here F = K, Bj, D) decays, the p parameter and the B — > X s 7 decay. In fllPf , all these 



processes have been analysed and two possible scenarios are obtained depending on the choice 
whether the constraint from R e b xp is enforced or not. Although the new experimental results 
are near the SM result, R% M = 0.2156, the world average value for R b (= 0.21656 ± 0.00074) is 
still almost one standard deviations higher than the SM one. This brings the possibility that 
an enhancement to the SM value is still necessary to get the correct experimental one. Such an 
enhancement is reached under the conditions ^ >> 1 and m H ± ~ 400 GeV |TtJ, where £^ is 
a model III parameter (see section 2) and v is the vacuum expectation value of the Higgs field 
responsible for the generation of fermion masses. 

First, the constraints for the FC couplings from AF = 2 processes for the model III were 



investigated withouth QCD corrections, under the following assumptions [K| 



1. Xij ~ A 



2. X U j = A$ « 1 , i,j = 1,2,3 , 

where u(d) is up (down) quark and i,j are the generation numbers. 

3. case 2 and further assumption 



X bb , X sb » 1 and X tt , X ct « 1 . (27) 



In the analysis, the ansatz 



«S d = Vt*, (28) 



is used. Note that this ansatz coincides with the one proposed by Cheng and Sher. 

Using the constraint coming from R e b xp , the measurement Br(B —>■ X s j), AF = 2 mixing 
and the result coming from the analysis of the p parameter, the following restrictions are 
obtained ||10|| : 

150 GeV < m H ± < 200 GeV , 
X bb » 1 , X tt « 1 , 

A 8b »l,A rf «l. (29) 
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Since the experimental results for R b p are still unclear, we disregard the constraint coming 
from R e b xp and we analyse the restrictions for the couplings ^ bb and ^ tc in the NLO 
aproximation, respecting the constraints due to the AF = 2 mixing, the p parameter and using 
the measurement by the CLEO M Collaboration: 



Br(B^X sl ) = (3.15 ± 0.35 ± 0.32) ■ 10 



-4 



Here, we explain why we use only the CLEO data in our analysis but not ALEPH one (Br{B — > 
X s 7) = (3.11 ± 0.80 ± 0.72) • 10~ 4 ). The ALEPH data has a larger error compared to CLEO 
data and it leads to a wide restriction region for \C% \, which includes the one coming from 
the CLEO data. Therefore, the CLEO data allows us to get more stringest constraints on the 
model parameters. 

The idea in this calculation is to take ^Ntc « ^mt^Nbb an d Cmb ~ QiCmj ~ 0, where 
the indices i,j denote d and s quarks. This choice permit us to neglect the neutral Higgs 
contributions because the Yukawa vertices are the combinations of £§ ib and ^ij- 
To reduce the b-quark mass dependence let us consider the ratio 

Br(B -> X sl ) 



H Br(B -> X c eu e ) 

\V; s V tb \ 2 Qa em //a> (3() , 



\V cb \ 2 ng(z)K(z 

where g(z) is the phase space factor in semileptonic b-decay, k(z) is the QCD correction to the 



semileptonic decay width [19 



g(z) = 1 - 8z 2 + 8z 6 - z 8 - 2AzHn z , 

= , _ 2o^) 2 _ 31 _ + 3 25 _ a _ 4 (l-^ ) 



and 2 = m c /m b . 



Using the CLEO data and following the same procedure as given in [13|], we reach the 
possible range for \C?\ as 

0.257 < \Cj ff \ < 0.439 . (32) 

In fig. ([TJ), we plot the parameter ^ tt with respect to £?} bb at n — 4.8 GeV and m H ± = 
400 GeV. We see, that there are two different restriction regions, where the upper one corre- 
sponds to the positive value, however the lower one to the negative value. Increasing 
^Nbb causes \£nu\ ^° decrease in both regions. With the given value of ^Nbb >> 1, the con- 
dition \r tb \ = < 1 is obtained. In the lower region it is possible that the ratio becomes 

10 



negative, i.e. r tb = < 0. Further, increasing m#± causes to increase |r t &| and the area of 

^N,bb 

the restriction region. 

Fig. (0) is devoted the same dependence as in fig. (HD and shows that the third region, which 
is almost a straight line, appears. In this region the ratio r tb » 1 and increases with increasing 
rriH± similar to the previous regions. 

Finally, we consider ^ tt dependence of ^ tc , which is a neutral FC coupling. In fig. (|3|) we 
plot the £jy u dependence of £j^ te for fixed £® bb = 60 m&, at fi = 4.8 GeV, and charged Higgs 
mass m H ± = 400 GeV. Here the selected region for ^ tt is 40 < ^% tt < 48. Increasing ^ tt 
forces the ratio r tc = 7^ to be negative. It is realized that the ratio \r tc \ becomes smaller 
when m H ± is larger. 

Still there is a region in which £jy te is constrained for the possible large value of ^ tt , namely 
£% >tt = 8.0 10 4 for m H ± = 400 GeV: 

- 0.24 <^ itc < 0.24, or 

-3.26 < £% ttc < -3.19 , (33) 

In conclusion, we find the constraints for the Yukawa couplings £^ it , £,§bb an d ^Ntc using 
the CLEO measurement Br(B — > X s 7) and respecting the restrictions due to the AF = 2 
mixing and the p parameter (see |10[ for details). The constraints for the other parameters of 



the model III from the existing experimental results require more detailed new analysis. 
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400 50-fi- — -*00 — — im — — 800 STVtt 



Figure 1: t^ tt as a function of ^ bb for the fixed value of the charged Higgs boson mass m^± = 
400 GeV at /i = 4.8 GeV. Here the constraint region is lying in between solid (dashed) curves. 
The solid (dashed) curves are the boundaries of the constraint region corresponding to > 

(c 7 e// <o) 
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Figure 2: Same as fig 1, but the third possible constraint region. 
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Figure 3: dependence of ^ ct for the fixed ^ bb = 60 mj, at // = 4.8 GeV and m fl ± = 
400 GeV. Here the constraint region is lying in between solid curves (dashed curves). The 
solid (dashed) curves are the boundaries of the constraint region corresponding to > 

(C 7 e// <0) 
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